Optimal quantum state reconstruction for cold trapped ions 
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We study the physical implementation of an optimal tomographic reconstruction scheme for the 
case of determining the state of a multi-qubit system, where trapped ions are used for defining 
qubits. The protocol is based on the use of mutually unbiased measurements and on the physical 
information described in H. Haffner et. al [Nature 438, 643-646 (2005)]. We introduce the concept 
of physical complexity for different types of unbiased measurements and analyze their generation in 
terms of one and two qubit gates for trapped ions. 
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A main task in any experimental physical setup for 
implementing quantum computation is the ability to de- 
termine the output state of any given quantum algorithm 
The standard procedure applied for quantum state 
reconstruction of a density operator lying in a 2^ dimen- 
sional quantum system, in the case of N qubits, consists 
in projecting the density operator onto 3^, completely 
factorized, bases in the corresponding Hilbert space 
All these measurements are obtained by applying rota- 
tions on single qubits (which are referred to as local op- 
erations) followed by projective measurements onto the 
logical basis. This was recently achieved for the case 
of eight qubits, with trapped ions The experiment 
was done by following the quantum computer architec- 
ture based on ions in a linear trap proposed by Cirac and 
Zoller Besides, the experimental implementations of 
several quantum protocols have also been reported by us- 
ing trapped ions [5] . In all these cases the quality of the 
protocols is tested using standard tomography for quan- 
tum state determination. This scheme has also been used 
in the cases of considering optical setups ^61] and NMR 
0. 

As it was mentioned above, in the standard measure- 
ment scheme only local operations are required to gen- 
erate all the necessary projections. In each basis (setup) 
2^ — 1 independent measurements can be performed, so 
that not all the experimental outcomes obtained in differ- 
ent bases are linearly independent, that is, there are re- 
dundant measurements. In the case of a A^-qubit system 
the anti-diagonal elements have the larger errors. Ac- 
tually, accumulated errors are not uniform; these errors 
depend on the number of single logic gates used for deter- 
mining given elements, so that larger errors appear when 
single logic gates act on all the particles. Assuming that 
there is an error e in the measurement of ion populations, 
then the accumulated error for anti-diagonal elements is 
of the order of £^2^-1-^2^-2(2^ - 1). These errors 
may lead to a density operator which does not satisfy 
the positiveness condition and so the information from 



the experimental data must be optimized. For this pur- 
pose the maximum likelihood estimation (MLE) method 
:8] has been used for the improvement of the density op- 
erators in experiments with light qubits as well as in 
experiments with matter qubits fs']. 

It is well known that the optimal quantum state de- 
termination is related to the concept of measurements 
on Mutually Unbiased Bases (MUB) which we will 
simply refer to as MUB-tomography. Such bases pos- 
sess the property of being maximally incompatible. This 
means that a state producing precise measurement re- 
sults in one set produces maximally random results in 
all the others. The set of mutually unbiased projectors 
given by 

= n=l,...,2^, a = l,..,2^-fl, (1) 

where 

Tr (P(")pif'^) - 5c.o.'5nn' + ^(1 - 60..') (2) 

and X]n=i = (-^ denotes the identity), defines a 
complete set of projection measurements, in the sense 
that the measured probabilities = Tr{Pn'^^ p) com- 
pletely determine the density operator of the system: 

a— 1 n— 1 

The number of MUBs for A^-qubits is (2^-1-1), which 
is essentially less than 3^. The use of MUBs can rep- 
resent a considerable reduction in the time needed for 
performing the full state determination. For instance, 
recently the reconstruction of a quantum state codified 
in the inner states of eight ions in a linear trap was re- 
ported [H . In this experiment the reconstruction process 
takes more than 10 hours, because of the measurement in 
6561 different bases and a hundred of times for each one, 
so the number of measurements is about 656,000, which 
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quickly grows for an increasing number of qubits. In the 
case of using MUB-tomography, the number of measure- 
ment bases is only 257 for determining all the elements 
of the density operator associated with this state. This 
could reduce the experimental time, roughly speaking, to 
25 minutes only, where we have considered a linear inter- 
polation. Recently, an alternative reconstruction scheme 
has proposed, which is based on using pyramidal states 
for single qubits H, 13] ■ The implementation of recon- 
struction using pyramidal states requires measurement of 
iV-particle correlations 

In the case of MUB-tomography each coefficient Pan 
in Eq. Q has an error associated with the measurement 
of only one projector Pn"^ Hence, the error in each 
coefficient is essentially determined by the ability of pro- 
jecting the system onto Pri""*- In practice, such measure- 
ments are implemented by projecting the system onto 
the logical basis after performing a set of unitary trans- 
formations. Such transformations, due to their nonlocal 
features, can be decomposed into a sequence of single 
and nonlocal gates, so that the error in this reconstruc- 
tion is mainly associated with the quality of these logic 
gates. Of course, the reconstructed density operator can 
be subjected to the MLE method. 

The main shortcut of the MUB-tomography scheme 
is the experimental implementation of MUB projectors, 
which is related to the fact that any set of MUBs contains 
non-factorizable bases. The measurements on such non- 
factorizable bases require application of non-local gate 
operations, which currently can not be performed with 
fidelity 1. Hence, the natural question for optimizing the 
experimental implementation can be casted as: Which is 
the set of MUBs that requires the minimum number of 
nonlocal operations? 

To approach this problem, we consider a complete set 
of MUBs where the basis factorization is denoted by 
the following set of natural numbers: {ki,k2, ...,fc0(jv)), 
where 4>{n) is the number of possible decompositions 
of 2^ + 1 as a sum of positive numbers, such that 
J2j^j =2^-1-1. Here, the notation {ki.k^,---) means 
that there are ki completely factorized bases; ^2 bases 
with two particle entanglement and all the other parti- 
cles are factorized, etc. Of course, for a given number of 
qubits, only a certain factorization structure is admitted. 
For instance, in the two-qubit case the only allowed struc- 
ture is (3,2), which means that there are 3 completely 
factorized bases and 2 non-factorized ones. In the case of 
3 qubits, the bases can be described by the following no- 
tation (n/, Hf,, n„/), with nf+n{,+nnf = 2^-1-1. Here, nf 
denotes the number of completely factorized bases; rib the 
number of bases with bipartite entanglement, i.e., each 
state of such bases is factorized as \ip)ij\ip)k,i j ^ k, 
where i,j,k — 1,2,3; and n„f represents the number 
of non factorized. For 3 qubits, there are 4 different 
sets of MUBs IJ, |l5|, which are denoted as (3,0,6), 



(1,6,2), (2,3,4) and (0,9,0). Because any 7V-qubit en- 
tanglement operation can be decomposed into a sequence 
of single qubit gates and controUed-NOT gates (CNOT 
gates) [13] , each basis (in the given set of MUBs) can be 
characterized by the minimum number CNOT gates. 

In this work we shall concentrate on the experimen- 
tal implementation of MUB-tomography in the case of 
trapped ions as reported in Ref. [5|]. In the case of 
'^°Ca^ trapped ions, where the qubit is codified in the 
ground state JO) = 5*1/2 (m = —1/2) and the metastable 
1) = D^/2{m = —1/2) state, single logic gates are im- 
plemented with a fidelity, ^sG: higher than 99%. How- 
ever, non-local operations are less accurate; here we will 
assume the reported value for the CNOT gate fidelity, 
'I'CNOT, which reaches a value up to 92(6)% [^. Then 
we can characterize the physical complexity of each set 
of MUBs as a function of the number of non-local gates 
needed for implementing the projection measurements. 
Then, a complexity of a given MUB is characterized by 
the number Cj, which can be defined as: 



C, = In ■ 



(4) 



CNOT 



where nj is the number of CNOT gates needed for gen- 
eration of such a basis. The total complexity of a given 
set of MUBs is then 



(5) 



Thus the total complexity is proportional to the total 
number of CNOT gates for preparing all the required 
projection measurements. As a simplest example let us 
consider the case of 3 qubits, where there are only 4 dif- 
ferent sets of MUBs. The bases labeled by (n/,n;,,n„/) 
have complexities C oc 0=i=n/-|-l*nb-l-2* ri„/. This 
means that the most adequate set for MUB-tomography 
is (0, 9, 0). We have assumed, in this derivation, that the 
fidelity of the factorized bases is one, in the sense that a 
sequence of single ion gates are required for their gener- 
ation. We also assume that CNOT gates between neigh- 
boring ions have the same fidelity as between ions that 
are further apart, because CNOT gates are implemented 
by using the center of mass motion as a data bus. The 
decomposition of the corresponding bases can be gener- 
ated by starting from the standard computational basis 
is that given in Table HI 

In the four-qubit case there are 34 sets of MUBs 
with different factorizations, which are labeled by 
(n/, rtf,, UtyUbb, rinf), with nj+nb+nt+Ubb+rinf = 
Then, there are Uf factorized basis; Ub basis with bi- 
partite entanglement, ij\ip) k\^) i] fit basis with tripar- 
tite entanglement, \il^)ijk\'p)i] nbb basis with two bipar- 
tite subsystems, \ip)ij\(p)ki; and n„/ non factorized basis, 
where i ^ j ^ k ^ I and i,j,k,l = 1,2,3,4. All such 
sets can be obtained in a regular way by applying finite 
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TABLE L Decomposition of the MUBs in the case of 3- 



qubits for (0,9,0) structure, where i?^-'' = Rl^'{n/2), with 

k = x,y,z are the single ions operations [l^ and Xcnot 
the controlled-NOT gate with the t-th and j-th ions as source 
and target, respectively. 



Basis 


Gate operations 


Basis 


Gate operations 


1 
2 
3 

4 
5 


p(l)p(3) (S3) p(2) 
iXx ACNOT^y 

p(l) n(3)^(23) p(2) 
J^y ACNOT^M 
p(3)^(23) n(2) 

p(3)^(13) p(l) 
^y XCNOT-'^a 
r(2)r(3) (13) „(1) 
J^x J^^x ACNOT 


6 
7 
8 
9 


r(2)r(3) (13) p(l) 
^y -f^z XCNOT-'^!/ 
p(2)p(3) (12) p(l) 
Jl-z -fix ACNOT^!/ 

XcNOT^!/ 
n(2)^(12) p(l) 
-fl-a: XCNOT-'^!/ 



phase-space methods [2J]. We remark that always there 
exists the so called "standard" set of MUBs, which is re- 
lated to rays in the finite phase space and can be easily 
constructed starting with two classes of operators con- 
taining either Oz or ax operators [2l|. The 4 qubit case 
is essentially different from the above discussed 3 qubit 
case. The main difference consists in that now there ex- 
ist two locally nonequivalent completely non-factorized 
states. Such states are isomorphic to the so-called graph 
states [17|. The first type. A, are isomorphic to the eigen- 
states of the set {azCTxII, OxOzOxl, Ictx^zCTx, Ilaxaz}- 
The second type,B, are isomorphic to the eigenstates of 
the set {az^xll, ax^z^x^x, IoxCTzI, IuxIctz}- Neverthe- 
less, both types of sets in the optimal decomposition are 
obtained by applying 3 CNOT gates only. Taking into 
account that MUBs with factorizations (1,3) and (2,2) 
can be generated with 2 CNOT gates, we realize that the 
optimum set of MUBs corresponds to the factorization 
structure (0, 0, 12, 2, 3), which only contains type A graph 
states. The complexity of such a set is C oc 37 while, for 
instance, the standard set of MUBs, (3, 0, 0, 2, 12), has a 
complexity C (x 40. 

The optimum set of MUBs corresponds to a set of 255 
disjoint operators which are arranged in a table consist- 
ing of 17 lines, so that each line contains 15 commuting 
operators 15|. The whole table can be obtained from 



only 8 elements arranged in two lines of commuting op- 
erators. 



la^^xC^y 


(Txaxlax 


l&zayCTz 


d'yazCTzax 


11a zl 


dy&xlax 


(^z(^z^zl 


axaxlaz 



(6) 



with factorization (4) and (1,3) for the first and second 
row, respectively. All other operators of the above ta- 
ble can be generated by using the following rule: Ar^c = 

* Ar^c-A and Ar^c = A2^c * Ai^r+c-3 for r > 2, 

where the index r (c) labels the r-th MUB (c-th opera- 
tor which has the basis states as eigenvectors of the r-th 
MUB), and r (c) goes from 1 to 2'* 1 (1 to 2'* - 1). 
Here, the sums are modulo 15. The four qubit case is 
somewhat special, because the standard table contains 
two bi-factorized bases (2, 2) and so there is no substan- 
tial difference with respect to the optimum set. This 



TABLE IL The decompositions on single logic gates and con- 
trolled NOT gates for the optimal set of MUBs (0, 0, 12, 2, 3) 
in case of 4-qubits. 



Basis 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 



Gate operations 



"* Xcnot Xcnot -^^ 

p(2)^(31) p(3)^{41) p(4) 

XCNOT^a; ACNOT^a; 
p(2) p(l)^(13) (14) p(l) 
J^v rLx XcnOtXcnOT-'^!/ 



R(4)p(2) (12) pii;^Aij; Ri 
n.x rty XCNOT^^: AcNOT^: 



(1),X13) 



p(l) p(2) (34) yyyO) \^-±) TJ 

Jl-t/ nx XcNOT^a; XCNOT^; 



(3), (24) 



p(4) p(3) (23) pi,^] 1,01; p 
n^y ^x XcNOT^a; AcNOT^: 



(2), (31) 



j(3) 



^y ' Xcnot Xcnot -'^^ 

r(3) r(1)^{14) p(l) (12) „(1) 
-fl-H ^y XCNOT-'^^; AcNOT-'^a 
r(1)r(2)p(4) (23) (24) p(2) 
J^x J^x J^x ACNOxXCNOT^y 
p(l)^(34) r(3)^{14) p(l) 
^x XCNOT^V XCNOT-'^a; 



(12) 



p(3) (24) r/(Z) yi-i) p^ 
XCNOT-^-H XCNOT^a; 
r(3),,(23) ,,(24) j,(2) 

ly rXy 



j(2), (14) 



j(l) 



-flu -flu Xcnot Xcnot -fly 



XCNOT^B XcNOT^a; 



-R 



E 



(2)^(14) p(l) (23) p(2) 

y XcNOT-^-a XcNOT^y 

(1)^(34) r(3) (32) (31) 
X XCNOT-^-a: XcNOtXcNOI 
(3),(13) p(l),(34) p(3) (12) 



T)K^)^K'-^) py'-i^y-''^) u^')^^'^'! rW 

^x AcNOT-'^a: XCNOT^^ XCNOT^B 
p(2)p(3) (12) p(l) (23) (24) p(2) 
J^x -fta: ACNOT-fiK XCNOTXCNOT-fl3= 



apparently does not occur for larger number of qubits, 
i.e, the standard table for TV > 4 would contain 2^ — 2 
completely non-factorized bases. The decompositions for 
the optimal set of MUBs is given in Table HIl 

The situation is quite different in the 5 qubit case, 
in which there exist at least 9000 non-isomorphic sets of 
MUBs with different factorizations, which are labelled by 
{nf,nb,nt,nbb,nf,nbt,nnf), where n/ denotes the num- 
ber of bases with four-particle entanglement, \ip)ijki \^)n, 
and Ubt the number of bases with bipartite and tripar- 
tite entanglement, \^)ijk\^)im, with i^j^k^l^m 
and i,j,k,l,m — 1,2,3,4,5. In case of 5 qubits there 
are four locally nonequivalent completely non-factorized 
states [ISj. In this case the standard table does not 
contain partially factorized bases and has the structure 
(3, 0, 0, 0, 0, 0, 30), and it is given by: 



llazCTzl 


^zf^z^zl^z 


llllaz 


l&zlazl 


lllazl 


ll&xaxl 


OxOxOx^^x 


llllOx 


l&xlaxl 


111-7x1 



The set of MUBs corresponds to 1023 disjoint operators 
which are arranged in a table consisting of 33 lines, so 
that each line contains 31 commuting operators. All 
the other projectors of this table are obtained by us- 
ing the following rule: Ar^c+is = Ar^c * ^r,c+i and 
Ar^c = Ai^c * A2^r+c-2, for r > 2, where the sums are 
modulo 31. This table contains 30 non-factorized bases 
among which there are three types corresponding to dif- 
ferent graphs: 6 of type B, 18 ofiype C and 6 of type 
D (graphs 6, 7 and 8 in Fig.4 [ll|). It results that 4 
CNOT gates are required to generate bases of type B 
and C, nevertheless the minimum number of CNOT gates 
to obtain the type D graphs is 5. This means that one 
needs 126 non-local operations to generate the whole set 
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of MUBs, with the corresponding complexity C oc 126. 
The optimum set has the structure (0,0,1,3,10,2,17) 
corresponding to a complexity C oc 112, and contains 
one non-factorized basis of type B and 16 bases of type 
C. In this case the set of operators needed for generating 
the whole table is 



ly x(y zfy zfy z 


axCyz^OyCTy 


(TzfTzllfTa; 


(y yO z<y z(y y(y z 


BxIfyzfyxOz 




OzOyOybzOx 


OzOz^OzOz 


G z ^GxGxO'x 


ioz (Ji^ii 



The rule for generating the other projectors is the same 
as in the standard table. For this table the decomposi- 
tions on single logic gates and CNOT gates can be ob- 
tained in the same way as for the above discussed cases. 

We have studied the problem of optimal tomographic 
reconstruction of a density operator of systems of N cold 
trapped ions. The optimality of a given set of MUBs is 
essentially defined in terms of the minimum number of 
required conditional operations. We have given an ex- 
plicit form of operations to generate the optimal set of 
MUBs for the case of three and four qubits, and the fac- 
torization operations for the five qubit case can be ob- 
tained in the same way. In the case of larger numbers 
of qubits, a generic procedure allowing generation of the 
whole set of MUBs is also available. It basically consists 
in finding all the possible strings of 2^ — 1 commuting 
operators (using explicit geometrical construction) with 
a subsequent separation onto the 2''^ -|- 1 disjoint sets [l^ . 
Eigenstates of such commuting operators in each disjoint 
set generate one of the corresponding MUBs [24] . This 
reconstruction scheme is valid for any physical setup, 
where non-local operations between neighboring qubits 
have the same fidelity as between distant qubits. This is 
satisfied in case of trapped ions because a CNOT gate, 
between any pair of ions, is implemented by using the 
center of mass motion as a data bus. If the physical im- 
plementation does not meet this requirement, the optimal 
MUBs to be used for reconstructing the state must have 
to be determined by considering the fidelities between 
non-neighboring qubits. 
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